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ABSTRACT 

With this Letter, we complete our model of the Galactic magnetic field (GMF), by using the WMAP7 
22 GHz total synchrotron intensity map and our earlier results to obtain a 13-parameter model of the 
Galactic random field, and to determine the strength of the striated random field. In combination 
with our 22-parameter description of the regular GMF, we obtain a very good fit to more than forty 
thousand extragalactic Faraday Rotation Measures (RMs) and the WMAP7 22 GHz polarized and 
total intensity synchrotron emission maps. The data calls for a striated component to the random 
field whose orientation is aligned with the regular field, having zero mean and rms strength ^20% 
larger than the regular field. A noteworthy feature of the new model is that the regular field has 
a significant out-of-plane component, which had not been considered earlier. The new GMF model 
gives a much better description of the totality of data than previous models in the literature. 



1. INTRODUCTION 

The magnetic field of the Galaxy (GMF) eludes our 
understanding in many respects, in spite of a great deal 
of effort measuring and modeling its various compo- 
nents. The most familiar components of the GMF are 
the large-scale regular fields a nd the smah -scale random 
fields, reviewed in Beuermann et al.|([l985|) and Sun et al. 
p008). The random fields are due to a number of phe- 
nomena including supernovae and other outflows, pos- 
sibly compounded by hydrodynamic turbulence, which 
are expected to result in randomly-oriented fields with a 
coherence leng th A of order 100 pc or less (Gaensler 



sity, for whi ch we adopt the standard NE2001 [Cordes fc 
Lazio (2002) model. The polarized and total synchrotron 



Johnston 1995| [Haverkorn et al. 2008) . In addition to 



these, there may be "striated" random fields whose ori 
entation is aligned over a large scale, but whose strength 
and sign varies on a small scale. Such striated fields can 
be produced by differential rotation of a medium contain- 
ing small scale random fields, the levitation of bubbles 
of hot plasma carrying trapped randomly oriented fields 
away from the disk, or both. Striated fields are a special 
case of the more generi c possibility o f aniso tropic ran- 
dom fields introduced in Sokoloff et al. (1998), which can 
be considered a superposition of multiple striated and 
pure l^y random fields. 



In Jansson & Farrar 



([2012) (JF12 below) we intro- 
duced a 22-parameter model of the large-scale, coherent 
or regular Galactic Magnetic Field, constrained by a si- 
multaneous fit to the WMAP7 Galactic Polarized Syn- 
chrotr on Emission (Stok es Q and [/, or, collectively, PI) 
maps (Gold et al.||201lD and more than forty thousand 
extragalactic I^araday rotation measures (RMs). The 
model includes an out-of-plane component and takes into 
account the possible presence of a striated random field 
and/or need for rescaling the assumed density of rela- 
tivist ic electrons. The new model gives a much better ac- 
counting of the observational data than achieved earlier. 
Key to the feasibility of the undertaking is the fact that 
the different data types probe different aspects of the field 
and are sensitive to different sources of uncertainty. The 
Faraday rotation measure is the integrated line-of-sight 
component of 5, weighted by the thermal electron den- 



mtensity are also line-of-sight integrals, but sensitive to 
the transverse rather than longitudinal component of the 
GMF and weighted by the relativistic (cosmic ray) elec- 
tron density, ncre- 

We present here a 13-parameter model of the spatial 
variation of the strength of the Galactic random field 
(GRF), obtained by applying the methodology and in- 
frastructure of JF12 to the total, unpolarized Galactic 
synchrotron emission from WMAP. The new random- 
field model provides a good fit to the large-scale structure 
of the total synchrotron emission / of the galaxy. This 
approach determines the rms value of 5rand as a func- 
tion of position within the Galaxy but is not sensitive to 
the "internal structure" of the random field, i.e., maxi- 
mum coherence length and power-spectrum. Combining 
these results with data on the variance can be used to 
constrain the distribution of coherence lengths, as will be 
presented elsewhere. 

A global parameterization of the GMF including both 
coherent and random fields is needed for many practical 
purposes. One important application is predicting the 
pattern of ultra-high energy cosmic ray (UHECR) de- 
fiections, to identify sources and gain insight into com- 
position; an early effort to characterize the Gala ctic ran- 
dom field for UH ECR studies was made by Tinyakov| 
& Tkachev (|2005|). Another important application is to 
understand the spatial variation and anisotropy of the 
diffusion coefficients for Galactic cosmic ray propagation, 
needed to obtain more realistic predictions for the diffuse 
gamma ray background and Galactic cosmic ray spectra 
and thus to have more accurate models for the back- 
grounds to possible Dark Matter signals. 

2. MODEL FIELD 

We allow for three distinct types of magnetic structures 
in our analysis: large-scale regular fields, striated fields, 
and small-scale random fields. These different structures 
can be disentangled because they contribute in general 
to different observables: the large-scale regular field can 
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contribute to all the observables, /, PI and RM; the stri- 
ated field can contribute to / and PI but not to RMs, in 
leading order, due to its changing sign; while the purely 
random field contributes only to the total synchrotron 
emission, /. 

In JF12 we considered two forms of ricre and adopted 
the GALPROP model provided by A. Strong (private 
communication, 2009), allowing for a possible overall 
rescaling by a factor a. Increasing the scale height 
and rescaling the GALPROP ncre, as could arise if 
anisotropic diffusion due to the out-of-plane field were 
incorporated i n GA LPR OP, further i mpro ves the fit; 

have also 



Strong et al.| ( |2011[ ) and |Jaffe et al. ( |2011| ) 



er i mpro ^ 
^ ( [2Qn| 
SHmcatior 



shown evidence ot a need tor a modification in the in- 
tersteller electron spectrum. In future work we intend 
to self-consistently solve for the GMF and ricre simul- 
taneously, but for the present we adopt JF12's partial- 
improvement of optimizing the normalization but not the 
shape of the GALPROP Ucre- 

The analysis of JF12 describes the coherent field as a 
superposition of disk, toroidal halo, and "X-field" com- 
ponents, with a total of 21 parameters. Each of these 
field components was initially allowed to have a separate 



amount of striation, with B'^^ ^ g^^- 



Pd^h^x Bl^^xreg and 



the orientation of a given component of the striated field 
taken to be aligned with corresponding component of the 
local coherent field. This introduces 3 more parameters 
than used to describe the coherent field alone. However 
there was not a significant improvement in the fit using 
the freedom of separate Pd,h,x^ so a single /3 value was 
used for all components in the final parameter optimiza- 
tion. When there is a single /3 for all components, it 
means that the possible striated field is aligned with the 
total local regular field resulting in a degeneracy between 
the strength of the striated component of the magnetic 
field and the relativistic electron density: if we write the 
multiplicative factor as 7 = a{l + we can interpret 
a as being a rescaling factor of the nominal relativistic 
electron density, with B^^^^ — PB^^^. The combined fit to 
Q, U, and RMs in JF12 established that 7 = 2.92 ±0.14; 
our present fit to / allows us to break the degeneracy and 
determine a and /3 separately. A future improvement will 
be to allow a and /3 to be spatially varying. 

After some experimentation, to get a good fit to the 
observations with a minimum of model complexity, we 
found that a good fit to the WMAP7 / map can be 
obtained by taking the random field to be a superpo- 
sition of a disk component, with a central region and 
eight spiral arms having the same geometry as for the 
coherent field, plus an extended random halo compo- 
nent. The total random field has an rms strength of 

-^rand 



B'^ 



The disk component of the GRF is modeled as the 
product of a radial factor times a vertical profile, taken 



to be a Gaussian of width Zq . In the inner 5 kpc, 
the radial factor is a constant: bint- Starting at 5 kpc, 
the rms strength is different from one arm to another. 
Its value at 5 kpc in the ith spiral arm is denoted 6^, 
and it decreases as ~ 1/r at larger radii. The extended 
random halo field is simply the product of an exponential 
in the radius and a Gaussian in the vertical direction, 
respectively: ^haio = ^0 exp[-r/ro] exp[-z'^ /2z^]. 




Figure 1. Top left: WMAP 22 GHz unpolarized synchrotron 
radiation. Top right: Simulated 22 GHz synchrotron radiation 
for the optimized model. Second row: As top row, in logarithmic 
scale. Third row, left panel: The estimated a of the unpolarized 
synchrotron data. Third row, right panel: The difference between 
observed and simulated data. Bottom left: The of the fit, with 
101 pixels removed (see section Isl. Bottom right: The of the 
masked-out pixels. Skymaps are m the Mollweide projection with 
galactic longitude zero in the center, increasing to the left. 



3. METHOD AND DATA 

Our approach to constraining the parameters of the 
GMF is described in detail in JF12 but we review it 
briefiy here, since the present analysis uses the same 
methodology. The observables used in the present anal- 
ysis are the a verage values of / in 13.4 square degree 
HEALpix cells dGorski et al.|2QQ5 ) and we minimize xj = 
X]i(^data,i " ^modei,i)^/c^/ i ^s a fuuctiou of the Galactic 
Random Field (GRF) parameters, taking the JF12 large 
scale GMF parameters and value of 7 as given. The con- 
ceptual foundation on which this and the JF12 project 
is based, is the recognition that the a^'s are themselves 
measurable observables. The variances in the primary 
observables are not merely due to observational uncer- 
tainties, but include and are dominated by the astro- 
physical variance caused by turbulent magnetic fields, 
inhomogeneities in the interstellar medium, and source 
contributions to the RMs. We measure (Ji for a given ob- 
servable, from the variance in values of that observable in 
the 16 sub-pixels of the zth 13.4 square degree HEALPix 
pixel. 

It is possible to estimate the unpolarized synchrotron 
emission without implementing a fully stochastic random 
field, which is computationally expensive, by generaliz- 
ing a trick introduced in Jansson & Farrar p012) for 
a striated random field. Assummg a spectraFTudex of 
3 for the relativistic electrons, the synchrotron emissiv- 
ity of a volume element in a region with local field B is 
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Figure 2. A histogram of for the 3072 pixels covering the 
sky. To reduce the impact of mainly nearby structures the 101 
pixels with x"^ greater than 10 are removed from the data, and a 
second optimization is performed to obtain the best-fit parameters 
presented in Table [l] 

(X Here B\ = Ticre ^^sin^^ where is the angle between 
the hne of sight and the direction oi B. If we add a ran- 
dom field of rms strength 5rand to ^^eg and average over 
the direction the random field may take, represented by 
< >, then < B\>= B'^^^s\v?e + | ^^and- To simulate 
the total synchrotron emission from a superposition of 
random, striated and regular fields, without ever creat- 
i n^ realizations of the random field, we use Hammurabi 
dWaelkens et aL|2QQ9D with the model 

ncre? but replacing 



reg, model 



1 



2 



)^reg,model 



(1) 

with a (1 + /3) =7 fixed from JF12. Since a and 1 + /3 
enter the expression for the total intensity differently, fit- 
ting the total intensity data allows them to be separately 
determined. As in JF12, we use a Metropolis Markov 
Chain Monte Carlo (MCMC) algorithm to find best-fit 
parameters and confidence levels. 

The WMAP total Galactic synchrotron intensity map 
used for the present analysis is shown in linear and log- 
arithmic scales in the top two rows of the left-hand col- 
umn in Fig. [l] and the map of measured a values for 
each 13.4 square degree HEALPix pixel is shown in the 
left panel of the third row. In order that fiuctuations due 
to nearby strong individual contributions do not distort 
the determination of the global properties of the Galac- 
tic random field, after performing an initial fit for the 
model parameters we mask-out the pixels with xf > 10, 
and then perform a second (and final) parameter opti- 
mization. The distribution of xf values from the initial 
model fit is shown in Fig. [2] The lower panels of Fig. [l] 
show the Xi per pixel maps (left) with the masked pixels 
removed and (right) showing the 101 masked-out pixels; 
the xf 's shown in these maps are from the second model 
optimization. 

4. RESULTS OF THE FIT 

The best-fit parameters for this model are given in Ta- 
ble [l] After finding the minimum chi-squared, an ad- 
ditional 100k MCMC iterations were performed to ob- 
tain a probability distribution in parameter space, re- 
fiecting the extent to which the observations do a better 
job constraining some parameters than others; the errors 
given in Table [l] are obtained from that distribution, by 
marginalizing over the other parameters. 



Table 1 

Best-fit parameters of the random field, with 1 - 



a intervals. 



Field 



Best fit Parameters 



Disk 

component 



b2 
b3 
b4 
b5 
be 
br 



^disk 



10.81 ± 2.33 /iG 
6.96 ± 1.58 /iG 
9.59±1.10/iG 
6.96±0.87/iG 
1.96 ± 1.32 /iG 

16.34 ± 2.53 /iG 
37.29 ± 2.39 /iG 

10.35 ± 4.43 /iG 
= 7.63 ± 1.39 /iG 



Description 



field strengths at r = 5 kpc 



field strength at r < 5 kpc 



Halo 


Bo 


= 4.68 ± 1.39 /iG 


field strength 


component 


ro 


= 10.97 ±3.80 kpc 


exponential scale length 






= 2.84 ± 1.30 kpc 


Gaussian scale height 


Striation 


^ = 


= 1.36 ±0.36 


striated field B'^^^. = (5B'^^^ 




Figure 3. Left panel: The random field in the disk. Right panel: 
The disk component of the JF12 coherent field model for compar- 
ison; it is clockwise in rings 3-6 and counterclockwise in 1,2,7, and 
8. 



The top two rows of Fig. [T] show the comparison be- 
tween data and predicted synchrotron intensity for this 
best-fit random field. The right-hand panel in the third 
row shows the residual between the WMAP 22- GHz to- 
tal synchrotron emission map and that produced by the 
random field model. It is concentrated near the Galac- 
tic plane and in isolated other regions, and is compatible 
with being due to individual sources such as Gen A, the 
North Polar Spur and other nearby structures not mod- 
eled here. The magnitude of the residual between model 
and data is small compared to the total signal, as most 
easily seen from the left panel of the second row, where 
the total intensity is shown in a logarithmic scale. 

The 13-parameter form of the Galactic random field 
adopted here proves to be sufficiently general to give a 
very good accounting of the data. In combination with 
the JF12 coherent field and the striated random field, it 
provides an excellent fit to the total intensity, /, with a 
reduced (x^ P^^" degree of freedom, Xdof ) ^^^^ ^^e sec- 
ond optimization of 1.064 with 2957 degrees of freedom. 
(One should not attach too much significance to the ex- 
act value of Xdof GRF model, since it depends on 
the arbitrary cut used to remove pixels with large "in- 
dividualistic" contributions; as evident from the map of 
residuals in Fig. [TJ additional pixels could be placed in 
this category which would decrease XdofO noted, the 
fit to / breaks the degeneracy between rescaling ncre and 
the presence of a striated random field. 

The left panel of Fig. |3] shows the disk component 
of the random field, with the magnitude of the coherent 
disk field from JF12 in the right panel for comparison. 
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The average rms strength of the disk component of the 
random field at the solar circle is 6.6 /iG, but it varies 
strongly from arm to arm. The spiral-arm model itself 
should not be taken too literally - it is presumably no 
more than a simplified encoding of some important fea- 
tures of the structure. 

Due to the large value of the random and striated fields 
compared to the coherent field, we cannot predict the 
value of the field at a particular position since the fiuctu- 
ations dominate the mean value. Nonetheless, we should 
check whether the predicted range of values for the total 
field in the solar neighborhood is consistent with observa- 
tions - keeping in mind that the rms component has 0(1) 
variance locally, so the actual local field at any particu- 
lar position can be expected to differ significantly from 
the result of combining the local coherent, striated and 
random components in quadrature. The estimate is addi- 
tionally uncertain due to our position near the boundary 
between the 4th and 5th arms, since the simple arm ge- 
ometry assumed in the present models is only expected 
to be valid in some average sense and may be locally 
modified, and the parameters specifying the geometry of 
the arms in JF12 were taken from the NE2001 model of 
Tie rather than being free parameters of the GMF model. 
With those caveats, the vertical component of the local 
coherent field is 0.2 /iG and the horizontal component is 
0.5-1.2 /iG in the 4th and 5th arms. Combining the co- 
herent, striated and random components in quadrature 
gives an estimate of the magnitude of the local field of 
3-5 /iG, with the range refiecting the values obtained for 
the 5th and 4th arms. Given the 0(1) variance from the 
random field, this estimated local GMF value is consis- 
tent with the 6 /iG value commonly ci ted (Beck 2008). 
The more re cent studies of [Taylor et al.| ( |20Q9p and |Mao| 
et al. (]2010l) are complementary, with the former hav- 
ing larger sky coverage and the latter a higher density 
of well-measured extragalactic sources but restricted to 
the polar caps. Our results are consistent with both, 
within the observationa l uncertaint i es and predicted fiuc- 
tuations. For instance, |Mao et al.| (|2010|) finds that the 
local random field is larger than the local coherent field, 
as we do, and estimates the halo random field above the 
solar system based on the variance in RMs in the po- 
lar caps to be 1 /iG in some average sense, consistent 
with our fit which decreases slowly from ^ 2 /iG in the 
Galactic plane, and is 1 /iG about 3 kpc above the plane. 

5. SUMMARY 

We have presented here a 13-parameter model of the 
random component of the Galactic magnetic field, to 
complete our characterization of the Galactic magnetic 
field. Taken together, our comprehensive GMF model fit 
utilizes 36 parameters, including the 21-parameter JF12 
model for the coherent field, the strength of a striated 
random field proportional to and aligned with the local 
coherent field, and an overall rescaling of the GALPROP 
cosmic ray electron density. This model provides an ex- 
cellent fit to more than 40,000 extragalactic RMs and the 
WMAP polarized and unpolarized Galactic synchrotron 
emission maps - a total of nearly 10,000 independent ob- 
servables, each with a measured variance. The striated 
random field has vanishing mean and is locally aligned 
with the regular field; its energy density is found to be 
^ 35% larger than that of the regular field. The fit in- 



dicates that the GALPROP cosmic ray electron density 
provided by A. Strong in 2009 is ^ 20% too low in the 
relevant energy regime, and/or its scale height may need 
to be revised; indeed, direct comparison to the current 
observations shows that including the correction factor 
inferred by our method gives a better fit in the relevant 
energy region. 

The GMF model provided here and in JF12 is a sub- 
stantial improvement over the previous state-of-the-art 
GMF models in several ways. First, the coherent com- 
ponent gives a much better fit to the RM, Q and U data 
than any other model. The JF12 coherent field model 
and 7 rescaling has a P^r degree of freedom of 1.096 
for the 6605 observables (pixels of RM, Q and U) while 
the fit is much worse for other recent m odels in the lit- 
eratur e, with Xdof 

of 2.663 & 4.971 for iPshi rkov et Z] 
(2011) BSS & ASS respectively and 1.672 lor Sun eTaT 
(2008), with their 2010 parameter update. Second, the 
functional form used in our model is better and more 
general: we enforce that the field is divergenceless and 
include out-of-plane and striated components. The im- 
provement due to our functional form is evident by re- 
optimizing the parameters of the other models with our 
MCMC procedure and Hammurabi; then their Xdof 
ues become 1.452, 1.591 and 1.325 respectively - still 
vastly worse than the 1.096 of JF12, considering there 
are nearly 6600 degrees of freedom. Finally, we pro- 
vide for the first time a detailed model of the random 
field strength. The final fit to 2971 pixels of / including 
the contributions of all three components of the GMF 
(coherent, striated and random) has Xdof ~ 1.064. We 
emphasize that there are two additional sources of uncer- 
tainty which we cannot quantify: the uncertainty in the 
relativistic (and thermal) electron distributions, needed 
to predict the RM and synchrotron emission observables 
for a given GMF model, and a possible inadequacy in 
the field model which could arise due to not having cho- 
sen a sufficiently general functional form. Future work 
will aim to address both of these issues, through self- 
consistent determination of the electron thermal and rel- 
ativistic electron distributions, and exploration of addi- 
tional field models. 

The total energy in the random and striated fields in 
this model is roughly six times that of the coherent field. 
Taking all three components together - random, stri- 
ated and regular - we find that the magnetic energy of 
the Galaxy is ~ 5 x 10^^ erg. The local magnetic field 
strength is predicted to be 3 — 5/iG, with an uncertainty 
of order of the rms local rando m field, 3/i G; this is con- 
sistent with the 6 /iG value of ( Beck||2008 ). The average 
rms strength of the disk component of the random field 
at the solar circle is 6.6 /iG, but has significant arm- 
to- arm variations. The vertical component of the JF12 
coherent local field is 0. 2 /iG, which is compatible with 
observational estimates ( Han fc Qiao||199^ Taylor et al. 
2009} [Ma^iriI][2MQl ) 



The improved determination of the Galactic random 
and coherent fields provided by this work represents a 
significant advance in our understanding of the Galaxy. 
Explaining the Galactic magnetic field will be an impor- 
tant challenge to theory, and using this more exact model 
will allow greater fidelity in making prediction for a host 
of important applications, from the defiections of ultra- 
high energy cosmic rays to the spectra and distribution 
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of Galactic cosmic rays. 
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